Fully fault tolerant quantum computation with non-deterministic gates 
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In certain approaches to quantum computing the operations between qubits are non-deterministic and 
likely to fail. For example, a distributed quantum processor would achieve scalability by networking 
together many small components; operations between components should assumed to be failure prone. In 
the logical limit of this architecture each component contains only one qubit. Here we derive thresholds for 
fault tolerant quantum computation under such extreme paradigms. We find that computation is supported 
for remarkably high failure rates (exceeding 90%) providing that failures are heralded, meanwhile the rate 
of unknown errors should not exceed 2 in 10 4 operations. 



The field of quantum information processing (QIP) has 
seen many experimental successes, but the challenge of 
scaling from a few qubits to large scale devices remains 
unsolved. One can argue that the issue is so crucial that 
it should dictate the choice of fundamental architecture for 
the machine. For example, in the concept of distributed 
QIP a plurality of small components, each similar in com- 
plexity to systems already realised experimentally, are net- 
worked together to constitute a full scale machine. The 
components may be trapped atoms or ions, or solid state 
nanostructures such as quantum dots or NV centres Q. 
Each component can be presumed to be under good con- 
trol, and it is understood that the key task is then to entangle 
the physically remote components. An attractive method of 
achieving this entangling operation (EO) is to arrange for 
each component to emit a photon that is correlated with 
the internal state of the component, before performing a 
joint measurement (with the aid of simple linear optical el- 
ements) of the photons. A considerable number of such 
entanglement schemes have been advanced since the first 
ideas in 1999 An important step was the realisation 

that photon loss can be detected, or heralded, within such 
a protocol |4][5). Generally in these remote entanglement 
protocols, one is supposed to employ optical measurements 
that simultaneously observe two, or even four 0, compo- 
nents simultaneously. This principle for generating entan- 
glement has in fact been demonstrated experimentally: first 
with ensemble systems and subsequently with individ- 
ual atoms |03- 

It is understood that the remote EOs may be failure 
prone. However, these failures are assumed to be heralded: 
the experimentalist is aware when a failure occurs. The ap- 
propriate strategy for dealing with such failures depends 
on the level of complexity within each component. In the 
case that each component incorporates multiple qubits then 
we can nominate one 'logical qubit' and use the other(s) to 
make repeated attempts at remote entanglement; when we 
are eventually successful then we can transfer the entan- 



glement to the logical qubits j9j[T0). However, many phys- 
ical systems may have only very limited complexity, and 
moreover it is always desirable to minimise the required 
complexity. Therefore it is interesting to consider the case 
of just one qubit in each component. This may be thought 
of as the extreme limit of the distributed paradigm. If we 
suppose that the probability of a heralded error is high, 
perhaps well above 50%, then it is clear that we cannot 
perform quantum computation by directly implementing 
a standard circuit model approach. However, it has been 
shown that even in spite of such heralded failures, arbi- 
trary quantum algorithms can be implemented J4j |5j [TD - 
[T5l . These insights are related to earlier ideas on photonic 
QIP |[l~6l[r71 . While such schemes demonstrated that large 
heralded failure rates can be tolerated, this was not shown 
in a fully fault-tolerant manner. In particular, it was not 
known if large heralded failure rates can be tolerated in the 
presence of realistic error rates for all other elementary op- 
erations. 

Fortunately, other studies have developed an approach 
which can be adapted to present purposes. Recently a se- 
ries of beautiful results by Raussendorf, Harrington and 
others described a method for QIP which involves creat- 
ing a large scale cluster state with a regular three dimen- 
sional lattice structure [ 1 8 - 20 ] . Defect regions within the 
3D lattice are braided together, yielding topologically pro- 
tected Clifford gates. QIP implemented using this topolog- 
ically protected cluster (TPC) state has a remarkably large 
tolerance against elementary errors (at rates < 1 %) during 
preparation, entangling operations and single qubit mea- 
surement. Subsequently, two of us have extended this idea 
to incorporate the possibility that the lattice contains a sig- 
nificant proportion of missing qubits at known locations 
(nearly 25% can be missing) ll2lTl23l . 

Here we consider the generation of a TPC state when the 
entangling operations are themselves subject to heralded 
failures during the cluster state growth process. The result 
is a lattice with a certain proportion of known failed entan- 
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glement relations (missing 'edges' in the graph state). The 
task of determining a threshold for universal QIP depends 
on proper choice of growth strategy together with a careful 
audit of the accumulation of unknown errors in that pro- 
cess. We show how to map this cluster state with missing 
'edges' to one with missing qubits, thereby making contact 
with the loss-tolerant thresholds quoted in the prior litera- 
ture EEHH- 

Several previous papers have considered the task of cre- 
ating large entangled states when the elementary EO is fail- 
ure prone (see Fig. [T] and caption). In principle a 'divide 
and conquer' approach can permit the entangled state to 
have positive growth on average for any nonzero success 
probability p s = 1 - p h Bl 151 [TTUT51 . Generally the so- 
lution involves generating relatively small resource states 
and subsequently connecting them. As shown in Fig. [TJa) 
the possible 'building block' resources include stars 11121 . 
linear clusters Bl fTTIl which in turn give rise to cross struc- 
tures |[T3lK and tree topologies ll24l . The last of these, 
also called the 'snowflake', has been proposed as a opti- 
mal choice for minimising errors 03). 

In the present paper our aim is to synthesise the TPC 
state (Fig. [TJ3, inset top left). This structure has the prop- 
erty that each node has four neighbours. Therefore we at- 
tempt to entangle together each resource with four others, 
as depicted in Fig. [lj>). In the particular example illus- 
trated we see that there will be N = 4- attempts to connect 
to each of the surrounding snowflakes. If one or more of 
these attempts succeeds, then have successfully connected 
the snowflakes, while with probability p% all attempts will 
fail and the resulting TPC state will have a missing 'edge' 
at that point. These missing edges are known, and there- 
fore are not errors but rather defects which we must allow 
for in the subsequent computation. Obviously, it will be 
necessary to create resource objects which are sufficiently 
large so that this net failure probability is below the thresh- 
old for fault tolerant QIP, which we presently discuss. For 
high values of ph we will see that the resource states must 
be considerably larger than those illustrated in Fig. [T] 

In order to evaluate this scheme, it is essential to deter- 
mine the accumulation of unknown errors when we per- 
form the star, cross, and snowflake strategies. These er- 
rors will occur during growth of the resource, during the 
fusion of resources (as in Fig.fTJ}), and also during the pro- 
cess of removing redundant qubits to simplify down to the 
TPC state. In order to minimise error accumulation during 
growth we make the aggressive choice that whenever there 
is a known failure during the growth of the resource object, 
that entire resource object is abandoned. Fortunately all 
three of the resources we consider - star, cross, snowflake 
- can be grown through a series of steps each of which 
(on success) doubles the entity's size. Thus the process is 
quick in the sense that it requires a number of successful 
steps that is merely a logarithmic function of the target re- 
source size. Here we assume each resource is grown using 
one of two forms of EO: we use either parity projections, 




FIG. 1: The figure shows graph states: nodes correspond to 
qubits, and connections ('edges') correspond to phase entangle- 
ment, (a) The 'building block' resources that have been consid- 
ered by previous authors, (b) Illustration of how one would syn- 
thesise the TPC state (depicted in upper left) by fusing together 
snowflake resource objects. One would use a quarter of the struc- 
ture to bond with each of four adjacent objects; basically the same 
approach applies for the star or the cross geometries. After the 
bonding stage, there must be a 'pruning' phase where we remove 
all but the core nodes and thus simplify down to the target lattice, 
(c) One minor revision in the case of the snowflake is that the fi- 
nal round involves fusing qubits that are one step from the 'core' 
in order to generate a 4-node. 



i.e. projecting a pair of qubits into the odd or even par- 
ity subspaces, or a canonical control-phase gate between 
the qubit pair, depending on which is more efficient J25j. 
Both operations are known to be possible through suitable 
measurements on emitted photons ll4ll5l [TD . 

We consider various forms of error. Single-qubit er- 
rors may occur during preparation, or while performing a 
single-qubit rotation, or during measurement. Moreover 
these errors may also occur passively in memory, i.e. there 
is a rate at which qubits decohere even when not part of any 
active operation. Meanwhile two-qubit errors may occur 
when we perform entanglement operations. We account for 
imperfections both in the emission of photons (e.g. from 
an imperfect selection rule in an atomic system, say) and 
errors arising from imperfect measurement of emitted pho- 
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tons ll25ll . Different probabilities are assigned to the various 
errors, however, for simplicity in generating the diagrams 
here we set the rates for all forms of the active 'gate' errors 
to be equal and we denote their probability pq. Memory 
errors are considered separately later in the paper. 

A principle conclusion from our analysis of error propa- 
gation is that two-qubit errors occurring during the growth 
and fusion of resource objects (e.g. snowflakes) typically 
appear as single qubit errors in the eventual TPC state. 
While there are instances where a two-qubit error can af- 
flict the ultimate lattice, the majority of these involve one 
qubit from the prime lattice and one from the dual, i.e. the 
black and green qubits in Fig. [TJ). Such correlations do 
not affect the fault tolerance threshold. There will be occa- 
sional instances of errors between two qubits both within 
the prime lattice, or both in the dual. However these are 
rare - for example in the case where one uses the snowflake 
strategy with pi, = 0.9, the rate for these errors is two orders 
of magnitude lower than the corresponding rate of single- 
qubit errors on the TPC (given equal rates for the various 
forms of error during growth) l25ll . In this case, two qubit 
errors only weakly affect the threshold lfl"9l. so as an ap- 
proximation we can ignore such events and assume that 
all gate errors affect at most one qubit in each sub lattice. 
Thus we consider a lattice with a (low) rate of random sin- 
gle qubit errors, and a (relatively high) portion of missing 
'edges' which are known. We need to determine the thresh- 
old for such a lattice to support computation. Fortunately, 
our previous work has considered the closely related case 
of a lattice with a significant number of missing nodes. We 
need only map the case of missing edges to that of missing 
nodes in order to make contact with that analysis and thus 
obtain thresholds in the present case. 

Consider the standard TPC state, specifically neighbour- 
ing qubits i (in the primal lattice) and j (in the dual lattice). 
Each qubit is centred on a face of its respective sublattice, 
and is a member of two cubic unit cells of the sublattice. 
In the ideal case where no bonds are missing, the product 
of cluster stabilisers associated with the faces of each cubic 
unit cell is simply the product of X operators acting on the 
respective face-centred qubits, yielding two parity-check 

1 2 

operators associated with each qubit: P i ' for qubit i and 

1 2 

P, ' for qubit j. Since these ideal parity check operators 
are just products of X operators on each face of the cor- 
responding cube, they commute point wise, which enables 
the error syndrome to be determined by single particle X 
measurements |fT8Tl2"0l . 

In the case where the bond between qubits i and j is 
missing, the cluster stabilisers associated with the missing 
bond are modified. Then, the product of cluster stabilisers 
centred on the cubic unit cell faces yields damaged parity 
check operators Pj' 2 = P}' 2 Zj and Pj' 2 = Pj' 2 Zi. Whilst 

A 1 2 "'12 

Pj ' and Pj ' commute, they do not commute point-wise 
(since [Xj,Zy] 7^ 0). In contrast to the ideal case, this means 
that determining the syndrome on the primal and dual lat- 



tices apparently requires measurement of the two-qubit op- 
erators XjZj and Z,X ; . 

Fortunately, by simply treating the qubits i and j at each 
end of the missing bond as though they were lost, and 
adopting the strategy in ETH%3"1 . we form products of the 
damaged parity check operators, yielding super-check op- 
erators Pi = P}Pf = PjPf and Pj = PjPj = PjPj. These 
new operators are independent of the qubits i and j, so they 
are unaffected by the missing bond between them. Fur- 
thermore each super-check operator involves only products 
of X operators from a single sublattice, so a missing bond 
manifests itself as a single missing qubit on each sublattice. 
This establishes a correspondence between missing bonds 
and correlated losses of neighbouring qubits. Error cor- 
rection is then realised by implementing the loss-tolerant, 
error-correcting protocol of 12TM231 to each sublattice in- 
dependently 11281 . 

Having made the connection to prior work on thresholds 
for the TPC state, we can now take any set of parameters 
for the low-level operations on qubits in the distributed ma- 
chine, compute the effective qubit loss rate, and determine 
whether quantum computation is possible. In Fig. [2] we 
show this phase diagram under the assumption that all gate 
error rates are equal. We see that very high rates of her- 
alded error can be tolerated, provided that the rate for un- 
known errors is below 2 x 10~ 4 . This is certainly a difficult 
number to achieve but might be possible in some imple- 
mentations, e.g. trapped ions for which multi-qubit mea- 
surements with fidelity around this rate have already been 
demonstrated 11261 . 

It remains to consider memory errors, which we assume 
happen at a lower rate than gate errors. In Fig. [3^ we show 
the effect of 'switching on' memory errors at a level equal 
to one tenth of the gate error rate. As one might expect, 
this lowers the overall threshold, but not dramatically. 

Finally we consider the question of physical resource 
scaling. From Fig.|2]one might be tempted to conclude that 
QIP is possible with extremely high rates of heralded error, 
perhaps reaching 99% or more. However, such a conclu- 
sion would neglect the ever increasing costs of preparing 
the resource objects. These objects become very large as 
Ph approaches unity. In Fig. [3]} we see that if pf, exceeds 
0.98, the size of each snowflake must be several thousand 
qubits. Recall that each snowflake ultimately corresponds 
to a single node in the TPC state, and therefore this factor 
would multiply the overhead already implicit in that ap- 
proach. However, values in the range of pi, w 0.9 may be 
tenable for technologies where the individual components 
of the distributed computer can be mass produced. 

In conclusion we have determined the threshold for 
quantum computation when two-qubit gates are non- 
deterministic. A specific case is that of a fully distributed 
machine, i.e. a network of components each of which 
contains only a single qubit. We find that it is tolerable 
if entanglement operations over the network fail with a 
rate exceeding 90%, provided that such failures are her- 
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FIG. 2: The principle results of our analysis. The lines define 
the parameter regimes where fault tolerant QIP is possible. Here 
for simplicity we set all gate errors, including both single-qubit 
and two-qubit errors, to be equally likely. This probability is de- 
noted pa and is plotted on the vertical axis. (Memory errors due 
to gradual decoherence are excluded, and are shown in Fig. 13) , 
Meanwhile, the probability pi, of an entanglement operation fail- 
ing in a heralded fashion is plotted on the horizontal axis. Note 
that p/ t can be very high, exceeding 90% if the snowflake strategy 
is employed. 



aided. The tolerable rate of un-heralded errors is 2 x 10~ 4 . 
Our analysis should allow experimentalists to determine 
if single-qubit components are feasible with their partic- 
ular approach, or if instead multi-qubit components must 
be adopted. 
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